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We study the electrophoretic mobility of spherical particles and the electrical conductivity in 
salt-free concentrated suspensions including finite ion size effects. An ideal salt-free suspension is 
composed of just charged colloidal particles and the added counterions that counterbalance their 
surface charge. In a very recent paper [Roa et at, Phys. Chem. Chem, Phys., 2011, 13, 3960- 
3968] we presented a model for the equilibrium electric double layer for this kind of suspensions 
considering the size of the counterions, and now we extend this work to analyze the response of the 
suspension under a static external electric field. The numerical results show the high importance of 
such corrections for moderate to high particle charges, especially when a region of closest approach 
of the counterions to the particle surface is considered. The present work sets the basis for further 
theoretical models with finite ion size corrections, concerning particularly the ac electrokinetics and 
rheology of such systems. 



I. INTRODUCTION 

The behavior of a suspension of charged particles un- 
der a static external electric field is a subject of study 
of electrokinetics and constitutes a classical problem 
H]. In 1978 O'Brien and White 3 revisited the prob- 
lem of electrophoresis and computed the electrophoretic 
mobility of a spherical particle for some cases of inter- 
est. A few years later O'Brien [Ij extended their work 
to obtain the electrical conductivity of a dilute suspen- 
sion of charged particles. Many of the classical stud- 
ies concern suspensions with low particle concentration, 
but nowadays the concentrated regime is the one that 
deserves more attention because of its practical applica- 
tions 01 . Ohshima developed analytic expressions for 
the electrophoretic mobility Q and the electrical con- 
ductivity [9| in concentrated suspensions by using a cell 
model approach. This approach has been successfully 
tested against experimental electrokinetic results in con- 
centrated suspensions From a theoretical point 
of view, these systems are difficult to understand due 
to the inherent complexity associated with the increas- 
ing particle-particle electrohydrodynamic interactions as 
particle concentration grows, and the possibility of over- 
lapping between adjacent double layers which will be un- 
avoidably present with high particle concentrations [ij] . 

In many typical cases, the presence of an external salt 
added to the system gives rise to an effective screening 
effect on repulsive electrostatic particle-particle interac- 
tions, depending on the salt concentration, which are 
mainly responsible, for example, for the generation of 
colloidal crystals or glasses. Thus, it would be of worth 
to study systems with a low screening regime for such in- 
teractions. Those systems are named salt-free because of 
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the absence of added external salt. The formation of col- 
loidal crystals is easier in this kind of systems, even at suf- 
ficiently low particle volume fractions |15l - [l8| . Of course, 
these salt-free systems contain ions in solution, the so- 
called "added counterions" stemming from the particles 
as they get charged, that counterbalance their surface 
charge preserving the clcctroneutrality. With the help of 
a cell model approach, Ohshima Il9i , and later Chiang 
et al. and Carrique et al. |2l| stressed the study 

of the electrophoretic mobility of spherical particles in 
salt- free suspensions. In the case of Carrique et al., the 
study was also extended to the computation of the elec- 
trical conductivity in salt-free concentrated suspensions. 
On the other hand, there is a lack of experimental results 
concerning the electrokinetic properties of salt-free sys- 
tems due to its difficult preparation. Between them, the 
experimental work of Palberg and coworkers jl6l - ll8l | is 
probably the most extensive using this kind of systems. 

All these theoretical studies are based on a mean-field 
description, that has a reasonable success when repre- 
senting the ionic concentration profiles at low to mod- 
erately charged interfaces. These studies also consider 
point-like ions, which, for highly charged particles, yields 
unphysical high counterion concentration profiles near 
such interfaces. In addition, these treatments neglect 
ion-ion correlations, which simplifies the real scenario. 
We can find in the literature different attempts to over- 
come these limitations. Some of them concern micro- 
scopic descriptions of ion-ion correlations and the finite 
size of the ions p3 - |26j that are able to predict impor- 
tant phenomena, like overcharging (g^j , but are basically 
restricted to equilibrium conditions. Others are based 
on macroscopic descriptions considering average interac- 
tions by mean-field approximations that include entropic 
contributions related to the excluded volume effect when 
the ions have a finite size [28l - [34| . The interested reader 
can find a historical overview about steric effects in the 
review of Bazant et al. jssj . Some simulations results 
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showed that these corrections work appreciably well with 
monovalent electrolytes for high surface charge densities 
and/or large ionic sizes [sHi- On the other hand, the 
macroscopic approaches permit us to make predictions 
under e quil ibrium and non-equilibrium conditions. Some 
authors [37h39I | have extended their works for equilibrium 
conditions to predict non-equilibrium properties, like the 
electrophorctic mobility or the electrical conductivity in 
diluted suspensions with electrolytes. 

One of the classical drawbacks of the mean-field ap- 
proaches without ion-ion correlations is their inability to 
explain important phenomena like overcharging, while 
Monte Carlo simulations achieve it by considering full 
ion-ion correlations. Very recently, Lopez-Garcia et al. 
[U, [3^ presented a modified standard electrokinetic 
model for diluted suspensions which takes into account 
the finite ion size and considers a minimum approach dis- 
tance of ions to the particle surface not necessarily equal 
to their effective radius in the bulk solution. They show 
that this model can predict overcharging in the case of 
high electrolyte concentrations and counterion valence. 
We think that this is a very important result because to 
our knowledge this is the first time that a phenomeno- 
logical theory based on macroscopic descriptions is able 
to predict this phenomenon. 

In a very recent paper [l^l we presented a model for 
the equilibrium electric double layer for spherical parti- 
cles in salt-free concentrated suspensions considering the 
size of the countcrions. The procedure was a general- 
ization of that already used by Borukhov [3l| for the 
special case of a salt-free suspension valid for the con- 
centrated case. Unlike Borukhov's treatment, our model 
also incorporates an excluded region in contact with the 
particle of a hydrated radius size, which has been shown 
by Aranda-Rascon et al. [s^ [33] to yield a more real- 
istic representation of the solid-liquid interface, and also 
to predict results in better agreement with experimental 
electrokinetic data. 

Our aim in this paper is to extend our previous work 
valid for equilibrium conditions to analyze the response 
of a salt-free concentrated suspension under a static ex- 
ternal electric field considering the size of the countcri- 
ons. We will study specially the electrophoretic mobil- 
ity of the particles and the electrical conductivity of the 
suspension. We will follow the treatment developed by 
Carriquc et al. plj for salt-free concentrated suspensions 
with point-like ions to achieve our electrokinetic model 
with ion size effects. 

The plan of this paper is as follows. In Section III Al we 
modify the governing electrokinetic equations to include 
the size of the countcrions and their distance of closest 
approach to the particle surface. The boundary condi- 
tions needed to solve the problem are discussed in Section 
IIIBI In Sections III CI and III Dl we present the expressions 
for the calculation of the electrophoretic mobility and 
the electrical conductivity, respectively. The results of 
the numerical calculations are shown in Section IIIII and 
analyzed upon changing particle volume fraction, parti- 




FIG. 1: Cell model including the distance of closest approach 
of the counterions to the particle surface. 

cle surface charge density, and size of the counterions. In 
order to show the realm of the finite ion size effect in salt- 
free suspensions, the results will be compared with the 
predictions that do not take into account a finite distance 
of closest approach to the particle surface, and also with 
the standard predictions for point-like ions. Conclusions 
are presented in Section IIVI 

II. THEORY 

A. Electrokinetic equations 

We use a cell model approach to account for the in- 
teractions between particles in concentrated suspensions 
through adequately chosen boundary conditions (bare 
Coulomb interactions among particles are included in an 
average sense, but ions-induced interactions between par- 
ticles as well as ion-ion correlations, are ignored). For 
details about the cell model approach see the excellent 
review of Zholkovskij et al. |4l| . In this approach, rep- 
resented in Fig. [1] each spherical particle of radius a is 
surrounded by a concentric shell of the liquid medium, 
having an outer radius b such that the particle/cell vol- 
ume ratio in the cell is equal to the particle volume frac- 
tion throughout the entire suspension, that is [i^ 



The basic assumption of the cell model is that the 
macroscopic properties of a suspension can be obtained 
from appropriate averages of local properties in a unique 
cell. 

Let us consider a spherical charged particle of radius 
a, surface charge density a and relative permittivity e^p 
immersed in a salt-free medium of relative permittivity 
ers and viscosity rj, with only the presence of the added 
counterions of valence Zc and drag coefficient Ac. In the 
presence of a static electric field E the particle moves 
with a uniform velocity Vg, the electrophoretic velocity. 
The axes of the spherical coordinate system (r, 9, if) are 
fixed at the center of the particle, with the polar axis 
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{9 = 0) parallel to the electric field. The solution of 
the problem requires the knowledge, at every point r of 
the system, of the electric potential, ^'(r), the number 
density of counterions, nc{r), their drift velocity, Vc(r), 
and the pressure, P{r). The fundamental electrokinetic 
equations connecting them are Q : the Poisson equation 
for the relationship between the electric potential and the 
charge density. 



V2vp(r) = -^ne(r) 



(2.2) 



the Navier-Stokes equation for low Reynolds number in 
the presence of an electrical body force for the fluid ve- 
locity, 

7/V2v(r) - VP(r) - Zcenc(r)V«'(r) = (2.3) 

the continuity equation for the counterions that implies 
the conservation of the number of counterions in the sys- 
tem, 



V- K(r)v,(r)] =0 



(2.4) 



and the Nernst-Planck equation for the flow of the coun- 
terions. 



nc(r)vc(r) ^ nc(r)v(r) - -^nc(r)V//c(r) (2.5) 

An 



where /ic(r) is the electrochemical potential of the coun- 
terions. We also take into account the continuity equa- 
tion for an incompressible fluid flow. 



V ■ v(r) = 



(2.6) 



In these equations, eo is the vacuum permittivity and 
e is the elementary electric charge. The drag coefficient 
Ac is related to the limiting ionic conductance or the 
difTusion coefficient Dc by 



AO 



Dr 



(2.7) 



where Na is Avogadro's number, fc^ is Boltzmann's con- 
stant, and T is the absolute temperature. 

As we are interested in studying the linear response of 
the system to an electric field, the following perturbation 
scheme is applied, where each quantity X is written as 
the sum of its equilibrium value, X", plus a perturbation 
term 5X linearly dependent with the field: 



*(r) = *°(r) +5*(r) 
71c (r) = ni{r) + Suciv) 
Aic(r) = ^2 + ^Mc(r) 



(2.8) 



As we showed in a previous paper [40[ , we introduce the 
finite size of the counterions by considering their excluded 



volume and including the entropy of the solvent molecules 
in the free energy of the suspension, F = U — TS 



U 



dr 



-|V*°(r) 



(2.9) 



TS = ksTnl 



dr 
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being 7i™°^ the maximum possible concentration of coun- 
terions due to the excluded volume effect, defined as 
j^maa: _ V^^, where V is the average volume occupied 
by an ion in the solution. The last term in eqn ()2.10|) 
is the one that accounts for the ion size effect, and was 
proposed earlier by Borukhov et al. [lOl ■ 

The variation of the free energy F = U — TS with 
respect to 4'°(r) provides the Poisson equation for the 
equilibrium 
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and the equilibrium counterions concentration is ob- 
tained performing the variation of the free energy with 
respect to n'^{r), obtaining 



,exp 



1 



cxp 



/ 2ce*0(r) A 
V kBT ) 



1 



(2.12) 



where be is an unknown coefficient that represents the 
ionic concentration where the electric potential is chosen 
to be zero. 

Wc also obtain the electrochemical potential doing the 
variation of the free energy with respect to rPci'^) [3^ . 
and assuming that the electrochemical potential out of 
equilibrium can be expressed in a similar way that in 
equilibrium conditions 



Aic(r) =2ce*(r) + fcsrin 




(2.13) 



In the case of equilibrium there is no external field 
and the particle is surrounded by a spherically symmet- 
rical charge distribution. Applying this symmetry and 
combining eqn (|2.1ip and (|2.12p . we obtain a modified 
Poisson-Boltzmann equation for the equilibrium electric 
potential 



d2^'0(r) ^ 2d*0(r) 
r 



dr^ 



dr 



ZcC- 



bc exp 



„e^°(r) 



eoCrs I _|_ 



exp 



- 1 



(2.14) 
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The electroneutrality of the cell implies that 



Q = 47ra^a = —Airzce i n°(r)r^dr 



(2.15) 



which is a necessary expression for the iterative cal- 
culation of the unknown he coefficient. An interested 
reader can find more details about the modified Poisson- 
Boltzmann equation for the equilibrium in ref. (40| . 

As indicated before, it is convenient to write the non- 
equilibrium quantities in terms of their equilibrium values 
plus a field-dependent perturbation. The symmetry of 
the problem allows us to define the functions h{r), 4'c{f), 
and Y{r) d 



v(r) ^ {Vr,Vg,V^) 



-hEcose, --^{rh)E sine, ] (2.16) 
r r ar 



(5/ic(r) ~ —Zce(f>c{r)E cos 6 



S'^i^r) = -r(r)^;cos^ 



(2.17) 



(2.18) 



with E = |E|. 

Substituting into the differential electrokinetic equa- 
tions, eqn (|2.2p - (|2.6p . the above mentioned perturbation 
scheme, neglecting nonlinear perturbations terms, and 
making use of the symmetry conditions of the problem 
we obtain 



C{Ch{r)) 



d*"(r) 
ksTrjr \ dr 

X n°(r) ( 0c(r) 



y(r) (2.19) 



the surface of the particle than their effective hydration 
radius, R, and, therefore, the ionic concentration will be 
zero in the region between the particle surface, r = a, 
and the spherical surface, r = a + R, defined by the 
counterion effective radius. This reasoning implies that 
counterions arc considered as spheres of radius R with a 
point charge at its center. 

With this consideration, the electrokinetic equations 
needed to solve the problem change into the following 
stepwise equations: the Poisson equation for the equilib- 
rium electric potential becomes 



2d^^0 a<r<a + R 

ar^ r ar — — 



(2.23) 



Eqn ([214]) 



R<r <b 



the Navier-Stokes equation for the fiuid velocity turns 
into 



CiCh{r))^0 a<r<a + R 



(2.24) 



Eqn (Pl^ a + R<r <b 



the equation for the conservation of the number of coun- 
terions now reads 



(j,^{r)=0 a<r<a + R 



(2.25) 



Eqn ^M) 



R<r <b 



and the Poisson equation for the perturbation of the elec- 
tric potential changes into 



CY{r) = a<r<a + R 
Eqn ([MT]) a + R<r <b 



(2.26) 



d*"(r) 



ksT \ dr 



Zc- 



d Mr) _ 2Ac h{r) 
dr e r 



CY{r) 



z^e^n°{r) 



(Mr) - Y{r)) 



where the C operator is defined by 

d2 2d 2 



C = 



dr-^ r dr 



(2.20) 



(2.21) 



(2.22) 



If we consider point-like counterions, n™°^ = oo, eqn 
(|2.14p . (|2.19p - (|2.21l) become the expressions obtained by 
Carrique et al. [2l|. 

Following the work of Aranda-Rascon et al. [s^, we 
incorporate a distance of closest approach of the coun- 
terions to the particle surface, resulting from their finite 
size. We assume that counterions cannot come closer to 



B. Boundary conditions 

We next specify the boundary conditions we use for 
the resolution of the electrokinetic equations. In the case 
of the equilibrium electric potential, we fix its origin at 
r — h, which results in 



*"(6) = 



(2.27) 



Using the electroneutrality condition of the cell, eqn 
(|2.15|) . and Gauss theorem to the outer surface of the 
cell, we obtain 



d^°{r) 



dr 



= 



(2.28) 



On the other hand, specifying the electrical state of 
the particle, and applying Gauss theorem to the outer 
side of the particle surface r = a we get 



d*°(r) 



dr 



(2.29) 
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We also force the equilibrium potential and its first 
derivative to be continuous at the surface r = a + R 
defined by the counterion effective radius. 

In the case of the electric potential out of equilibrium, 
the discontinuity of the normal component of the dis- 
placement vector at the particle surface of charge density 
a states 



,V*(r) 



—a 



(2.30) 



where \l/p(r) is the electric potential in the interior region 
of the solid particle, and f is the normal unit vector out- 
ward to the surface. Also, the continuity of the electric 
potential at the surface of the particle has to be consid- 
ered 



^ ' I r—a ^ ^ ^ \ r—a 



(2.31) 



According to Shilov-Zharkikh-Borkovskaya boundary 
conditions [4j] , the connection between the macroscopic 
experimentally measured electric field (E) and local elec- 
tric properties is 



*(r)| ,-*"(r)| 



r—b 



(E) 



I r—b 



(2.32) 



We also must impose the continuity of the electric po- 
tential out of equilibrium and of its first derivative at the 
boundary surface r = a + R. 

Following again Shilov-Zharkikh-Borkovskaya bound- 
ary conditions, the ionic perturbation at the outer surface 
of the cell must be zero, or equivalently 



"c(r)| 



(2.33) 



As the solid particles are impenetrable objects for the 
ions, the velocity of the ions in the normal direction to 
the particle surface is zero 



= (r) 



= 



(2.34) 



Due to the inclusion of a distance of closest approach 
of the counterions to the particle surface, the density of 
counterions nc(r) and their drift velocity Vc(r) will be 
discontinuous at the surface r = a + i?, being zero in the 
region [a, a + R] and non-zero in the region [a + R,b]. 

The liquid located at the particle surface, r = a, is 
considered immobile, strongly attached to the particle. 
The latter fact that the liquid cannot slip on the particle 
is expressed as 



v(r)| =0 

^ / I r—a 



(2.35) 



At the outer surface of the cell, r = b, we follow 
Kuwabara's boundary conditions [4J|. In the radial di- 
rection, the velocity of the liquid far from the particle 
will be the negative of the radial component of the elec- 
trophoretic velocity 



v(r) 



.(r) 



(2.36) 



According also to Kuwabara, the fluid flow is free of 
vorticity at the outer surface of the cell 



V X v(r) 



r=b 







(2.37) 



At the boundary surface r = a + i? we must consider 
the continuity of the normal and tangential components 
of the fluid velocity as well as the continuity of vorticity 



and pressure [45 



Finally, in the stationary state, the net force acting on 
the particle or the unit cell must be zero. Since the net 
electric charge within the unit cell is zero, there is no 
net electric force acting on the unit cell, and we need to 
consider only the hydrodynamic force. For details about 
this boundary condition see ref. Q or Appendix 1 in ref. 

M- 

In terms of the radial functions ^'"(r), Y{r), (f>c{r) and 
/i(r), the previous boundary conditions change into: 
(i) at the particle surface r = a 



d*°(r) 



dr 



dY{r) 



dr 



Yia) 



= 



h{a) = 



dh{r) 



dr 



(2.38) 



(2.39) 



(2.40) 



(2.41) 



(ii) at the surface r = a + R defined by the counterion 
effective radius 



*"(a + i?-) = ■i'°{a + R+) 



d*0(r) 



dr 



d*0(r) 



r=a+R- 



dr 



r=a+R+ 



Y{a + R-) = Y{a + R+) 



dY{r) 



dr 



dF(r) 



r=a+R- 



dr 



dr 



r=a+R+ 



r=a+R+ 



h{a + R^) = h{a + R-^ 



dhjr) 
dr 



r=a+R- 



dfe(r) 
dr 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



(2.47) 



{2.4 



r=a+R+ 
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Ch{a + R-) = Ch{a + R+) 



(2.49) 



d^h{r) 



r=a+R- 



dr^ 



r=a+R+ 



^^^^^^nl{a + R+)Y{a + R+) (2.50) 



{in) and finally, at the outer surface of the cell r = b 



^-"(6) = 



d*0(r) 



dr 



r—b 



Y{b) = b 



(2.51) 



(2.52) 



(2.53) 



The macroscopic electric field (E) is given by 
1 



(E) 



Vr 



cell J V^^ii 



V*(r)dF 



(2.60) 



The electric current density of a salt-free suspension is 
defined by 

J(r) = Zcenc(r)vc(r) 



= z,en, l^v(r) - — V/^c(r) j (2.61) 

Following a similar procedure to that described for the 
conductivity of suspensions in salt solutions in ref. [47j 
we obtain (see also ref. plj ) 



K = 



zle^ d(/>c(r) 



\r dr 



r—b 



Mb) = b 



(2.54) 



III. RESULTS AND DISCUSSION 



Ch{b) = 



(2.55) 



7?— [rCh{r)] - z,eb,Y{b) = (2.56) 



C. Electrophoretic mobility 

The electrophoretic mobility fi oi a spherical particle in 
a concentrated colloidal suspension can be defined from 
the relation between the electrophoretic velocity of the 
particle Vg and the macroscopic electric field (E). From 
the boundary condition, eqn p.36p . the definition |ve| = 
/i|(E)|, and the symmetry cqn (j2.16p . we obtain 



2h{b) 



As usual, the mobility data will be scaled as 

3r/e 



(2.57) 



(2.58) 



2eo€rskBT 

where n* is the nondimensional electrophoretic mobility. 

D. Electrical conductivity 

The electrical conductivity, K, of the suspension, is 
usually defined in terms of the volume averages of the 
local electric current density and electric field in a cell 
representing the whole suspension. 

(J) = / J(r)dy = Km (2.59) 



We will discuss the results obtained from three differ- 
ent electrokinetic models: the classical with point-like 
countcrions (PL) from ref. [2l|, its modification to take 
into account the finite ion size (FIS), eqn (j2.14p . (j2.19|) - 
(|2.2ip . and the complete electrokinetic model that also 
considers the distance of closest approach of the countc- 
rions to the charged particle surface (FIS-I-L), eqn (|2.23p - 
()2.26p . In this complete electrokinetic model the electric 
potential fulfills the Laplace equation in the excluded re- 
gion in contact with the particle. 

The governing electrokinetic equations with their 
boundary conditions form a boundary value problem that 
can be solved numerically using the MATLAB routine 
bvp4c [i^. This routine computes the solution with a 
finite difference method by the three-stage Lobatto IIIA 
formula, which is a collocation method that provides a 
solution that is fourth order uniformly accurate at all 
the mesh points. The resulting mesh is non-uniformly 
spaced out and has been chosen to fulfill the admitted 
error tolerance (always taken lower than 10~^). 

For all the calculations, the temperature T has been 
chosen equal to 298.15 K, the viscosity of the solution 
7] = 0.89-10~^ Pa-s, and the relative electric permittiv- 
ity of the suspending liquid 6,.^ = 78.55, which coincides 
with that of the deionized water, although no additional 
ions different to those stemming from the particles have 
been considered in the present model. We have used the 
value erp = 2 for the relative permittivity of the particles. 
Also, the particle radius a has been taken equal to 100 
nm and the valence of the added counterions Zc = -1-1. 
Other values for Zc could have been chosen. The model 
for point-like ions is able to work with any value of Zc, 
but we think that the predictions of this model will be 
less accurate in the case of multivalent counterions, since 
it is based on a mean-field approach that does not con- 
sider ion-ion correlations. Nevertheless, when we take 
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FIG. 2: Scaled electrophoretic mobility against the parti- 
cle surface charge density for different ion sizes, considering 
(dashed lines) or not (solid lines) the excluded region in con- 
tact with the particle. Black lines show the results for point- 
like ions. 



FIG. 3; Scaled electrophoretic mobility against the particle 
volume fraction for different ion sizes, considering (dashed 
lines) or not (solid lines) the excluded region in contact with 
the particle. Black lines show the results for point-like ions. 



into account the finite size of the ions, the main objec- 
tive of this work, we include correlations associated with 
the ionic excluded volume, solving partly this problem, 
because we are still not considering the electrostatic ion- 
ion correlations. 

For the sake of simplicity, we assume that the average 
volume occupied by a counterion is y = {2R)^, being 2R 
the counterion effective diameter. With this considera- 
tion, the maximum possible concentration of counterions 
due to the excluded volume effect is n^"^ = {2R)^^. 
This corresponds to a simple cubic package (52% pack- 
ing). In molar concentrations, the values used in the cal- 
culations, n™"^ = 22, 4 and 1.7 M, correspond approxi- 
mately to counterion effective diameters of 2R = 0.425, 
0.75 and 1 nm, respectively. These are typical hydratcd 
ionic radii [i^. The diffusion coefficient of the counteri- 
ons has been chosen as Dc = 9.34- 10~^ ni^/s, that cor- 
responds to the value for H"*" ions, which are commonly 
found in many experimental conditions with pure salt- 
free suspensions, although other different values could 
have been used. 



A. Electrophoretic mobility 

The classical behavior of the electrophoretic mobility 
of spherical particles in salt-free concentrated suspen- 
sions when we consider point-like counterions, PL model, 
is as follows: for low particle surface charges, there is a 
large increment of the electrophoretic mobility with the 
surface charge. When the particle volume fraction de- 
creases, this increment is even larger. This behavior sat- 
isfies a Hiickel law linearly connecting both magnitudes. 
When the particle charge increases, the electrophoretic 
mobility reaches a plateau and becomes practically inde- 



pendent of particle charge. This fact has been associated 
with the generation of a condensation layer of counterions 
close to the particle surface [l^ . Between these regimes 
there is a maximum followed by a small diminution of the 
electrophoretic mobility that depends on particle volume 
fraction. Also, the lower the particle volume fraction, 
the higher the electrophoretic mobility for every parti- 
cle charge value. These classical behaviors arc shown in 
solid black lines in Figs. [2] and [31 The results displayed 
in solid colored lines in Figs. |4] and [5] are also useful to 
understand this discussion. 

If we take into account the size of the counterions, 
FIS model, we find deviations from the point-like case. 
As wc can sec from the results of the solid colored lines 
in Fig. [21 the small diminution passed the maximum 
in the electrophoretic mobility tends to disappear when 
the ion size becomes important. Moreover, if the size of 
the counterions is sufficiently large, also the maximum 
disappears. In this case we find two different regimes for 
the electrophoretic mobility upon changing the surface 
charge of the particles: the initial large increment of the 
mobility with the surface charge, similar to the one for 
point-like ions, now followed, for higher particle charges, 
by another region with a small rate of increment of the 
electrophoretic mobility. 

When we study the behavior of the electrophoretic mo- 
bility when changing the particle volume fraction, we ob- 
serve that the results of the FIS model differ from those 
for point-like ions, obtaining higher values for the mobil- 
ity as we approach to the concentrated regime, see solid 
colored lines in Fig. [31 Also, if the ion size is sufficiently 
large, we can find a broad minimum at high particle vol- 
ume fraction, in contrast with the PL case. 

If the counterion size approaches to zero, or equiva- 
lently n™"^ oo, the results of the FIS model approxi- 
mate to those of the PL model in any case. Also, for low 
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-20 

a (|iC/cm^) 

FIG. 4: Scaled electrophoretic mobility against the surface 
charge density for different particle volume fraction values. 
Solid lines show the results for point-like ions. Dashed lines 
show the resuhs of the FIS-I-L model with n^""^ = 4 M. 




FIG. 5; Scaled electrophoretic mobility against the particle 
volume fraction for different surface charge densities. Solid 
lines show the results for point-like ions. Dashed lines show 
the results of the FIS-^L model with n'""'' = 4 M. 



particle charges and low particle volume fractions the re- 
sults of both models are nearly the same. For the remain- 
ing situations we always observe that the FIS model pre- 
dicts higher values of the electrophoretic mobility than 
those calculated with the classical model for point-like 
ions, whatever the ion size. 

If we also consider the distance of closest approach of 
the counterions to the particle surface, FIS-I-L model, we 
obtain large differences in comparison with the PL model 
from moderate to high particle charges, see dashed lines 
in Figs. [2]and|4l Although it is not shown in Fig. [2l when 
the particle surface charge is extremely high, the elec- 
trophoretic mobility predicted with the FIS and FIS+L 
models must reach the same plateau value, because the 




FIG. 6: Scaled polar component of the fluid velocity over the 
particle equator, 9 — 7r/2, (a) and scaled radial component of 
the fluid velocity at the front of the particle, 6 = n, (b) along 
the cell for different ion sizes. Different lines have the same 
meaning than those shown in Fig. [21 



distance of closest approach becomes negligible in com- 
parison with the width of the condensation counterions 
layer. 

When we change the particle volume fraction, dashed 
lines in Fig. [31 the consideration of the excluded region 
in contact with the particle augments the effect that we 
observed with the FIS model. The results displayed in 
Fig. [5l also show how for a fixed size of the counterion, 
the FIS-l-L model predicts large deviations from the PL 
model for concentrated suspensions (dashed lines versus 
solid lines). We observe that always the FIS-I-L model 
predicts equal or higher values of the electrophoretic mo- 
bility than the PL and the FIS models, whatever the ion 
size. 

As a main conclusion from Figs. [Utoini we can affirm 
that the consideration of finite ion size effects leads to an 
increase of the electrophoretic mobility over the PL case, 
using H+ ions as counterions. As the ionic concentration 
in the cell has been altered, we will analyze the changes 
in the convcctive fluid flow, the counterions fluxes and 
the perturbed counterions concentration, as well as the 
overall forces and polarizations induced by the electric 
field for the FIS and FIS-I-L models in comparison with 
the PL case. 

Fig. [SJi shows the scaled polar component of the fluid 
velocity over the particle equator, 9 = 7r/2, and Fig. |6)d 
presents the scaled radial component of the fluid veloc- 
ity at the front of the particle, 9 = tt, along the cell. 
Solid black lines, solid colored lines and dashed lines. 
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FIG. 7: Scaled polar component of the flux of counterions over 
the particle equator, 6 — it/2, (a) and scaled radial component 
of the flux of counterions at the front of the particle, 9 — tt, 
(b) along the cell for different ion sizes. Different lines have 
the same meaning than those shown in Fig. [5] In the inset of 
Fig. [7^ we enlarge the region close to the particle surface. 



represent the results of the PL, FIS and FIS+L models, 
respectively. Different colors stand for different counte- 
rion sizes. The particle surface charge density has been 
chosen equal to —40 /iC/cm^, and the particle volume 
fraction is (/) = 0.5, which implies a normalized cell size 
of b/a = 1.26. We define the scaled fluid velocity as 



2eQer sksTE 



v(r) 



(3.1) 



where the different components of v(r) can be obtained 
from eqn (|2.16p . Both quantities, vg at 7r/2 and Vr at 
TT, are of interest because they give us an idea of the 
magnitude of the clcctrophorctic velocity because they 
are antiparallel to it. 

We can see how the polar fluid velocity increases with 
the distance to the particle surface, with a high rate close 
to the particle surface, and diminishing it as we approx- 
imate to the outer surface of the cell. In the case of the 
radial component, there is a linear increase after an initial 
slower growth rate very close to the particle surface. This 
behavior is the same for the PL, FIS and FIS+L models. 
We also observe that the numerical values obtained with 
the FIS+L model are larger than those obtained with the 
FIS model, being the latter results also larger than those 
obtained with the PL model, in concordance with our 
predictions for the mobility. 

Fig. [7^ shows the scaled polar component of the coun- 
terions flux over the particle equator, 9 = 7r/2, and Fig. 
[7}d presents the scaled radial component of the flux of the 



counterions at the front of the particle, 6 = tt, along the 
cell. Different lines have the same meaning than in Fig. 
El We define the scaled counterions flux as 

J:(r)-x^^^^f^Jc(r) (3.2) 

being Jc(r) = nc(r)vc(r), where Vc(r) is obtained from 
eqn (|2.5p . We can see in Fig. [7] how the inclusion of 
finite ion size effects largely decreases the magnitude of 
the counterions fiuxes close to the particle surface in both 
the polar and the radial directions in comparison with the 
PL case. In addition, the counterions fiuxes are highly 
increased in the FIS and FIS+L cases as we move away 
from the particle surface because now both the coun- 
terions concentration and the counterions velocity (not 
shown for brevity) reach higher values. 

The enhancement of the fiuid velocity observed in the 
FIS and FIS+L cases in comparison with the PL model 
could be associated with the increment of the counterions 
fiuxes in the region not in the immediate vicinity of the 
particle, because the ions have been expelled out due to 
the excluded volume effect, and consequently the electri- 
cal body force in that region is greater. However, as we 
can see in Table HI the total electrical body force, —F*, 
is lower in the FIS and FIS+L cases, in comparison with 
the PL case. This is due to a very large diminution of the 
electrical body force very close to the particle when the fi- 
nite ion size is taken into account, not only because ions 
have moved away from the vicinity of the particle sur- 
face but also for the remarkable diminution of the local 
electric field due to the increased induced electric polar- 
ization (see Table H]). 

The total hydrodynamic and electric forces acting on 
the particle shown in Table U can be calculated by 



a" ■ f dS^ 



~M 



■ f dSr, 



(3.3) 



(3.4) 



TABLE I: Scaled hydrodynamic, F^, and electric, F* , forces 
acting on the particle in the direction of the field. The to- 
tal body force in the fluid is equal to ~F* due to the elec- 
troneutrality of the cell. (Pc)* and (Pd)* are the scaled in- 
duced charge and dielectric polarizations in the direction of 
the field, respectively. All the calculations were performed at 
a = -40/iC/cm^ and = 0.5. 



Model / n^"" 


f: 


f: 


{Per 


{Pny 


PL 


685 


-696 


346 


-71.7 


FIS / 22 M 


678 


-677 


350 


-71.8 


FIS / 4 M 


582 


-582 


368 


-72.6 


FIS / 1.7 M 


492 


-492 


399 


-73.4 


FIS+L / 22 M 


524 


-524 


365 


-73.1 


FIS+L / 4 M 


447 


-447 


386 


-73.8 


FIS+L / 1.7 M 


394 


-394 


418 


-74.2 
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where Sp is the surfaee of the solid partiele, and and 
CT*^ are the hydrodynamic and the Maxwell stress ten- 
sors, respectively Evaluating these expressions, we 
obtain 



Zce'nP^{a)Y {a) 



4 . Y{a) 
3 a 



^ + 2 

67- t- 



k (3.5) 



(3.6) 



where k points to the direction of the macroscopic electric 
field. Both forces are scaled as follows 



3e 



inaeoerskBTE 



(3.7) 



The numerical values of the F* are negative because 
this force has the opposite direction of the electric field for 
a negative particle. The hydrodynamic force opposes the 
movement of the particle and therefore has the direction 
of the field. In the stationary state the total force acting 
on the particle must be zero, as we can see by summing 
the values of both forces in Table HI The small numerical 
discrepancies observed, mainly for the PL case, can be 
removed by improving the mesh for the resolution of the 
electrokinetic equations very close to the particle surface 
although a large computational time is required. 

It is worthwhile to mention that both the hydrody- 
namic and the electric forces calculated in Table U for the 
FIS and FIS-I-L models are lower than those of the PL 
case. The electric force acting on the particle has a driv- 
ing and a relaxation contribution. As the driving force is 
constant in the study displayed in TablelH the diminution 
of the total electric force is related with a change of the 
relaxation force as finite ion size is taken into account. 

This relaxation force will depend on the electric dipole 
moment induced on the particle and its double layer by 
the electric field. A related quantity is the average in- 
duced dipole moment density whose components are: the 
charge polarization, (Pc), and the induced dipole mo- 
ment density arising from the polarization of the dielec- 
tric continuum of the medium and the particles, (Pd) 



{Pd) = 
where 



Pc) - (t^ I vz,e5n,{v)Av\ (3.8) 

\ Vcell Jv^^n I 

-1/ r(e(r) - eo)VJvI/(r)dA (3.9) 



e r 



(3.10) 



and Vp and are the particle and solution volumes in the 
cell, respectively. According to the procedure developed 
for AC electric fields and point-like ions by Bradshaw- 
Hajek et al. (soj . particularized to its DC limit and also 
accounting for the distance of closest approach of the 
counterions to the particle surface, the latter equations 
become 



(Pc) 




dV 
dr 



r=a+R+ 




(Cr 



+ (frs 



E) (3.11) 



(E) (3.12) 



We show in Table |T] the scaled average polarization 
contributions, (Pc)* and (Pd)*, calculated by 



(Pc.g) 
eoE 



(3.13) 



The charge polarization takes positive values and 
therefore the induced dipole moment generates an elec- 
tric field which opposes the external one, thus penalizing 
the particle movement. On the contrary, the induced 
dipole moment due to the dielectric polarization has the 
opposite direction, reinforcing the effect of the external 
electric field on the particle movement. As the charge 
polarization contribution is considerably larger than the 
dielectric polarization one, see Table HI the effect on the 
particle movement will be a net relaxation force that op- 
poses the electric driving force, yielding a smaller total 
electric force in the FIS and FIS-I-L models, in compari- 
son with the PL case (see Table |I| . 

In Fig. |8] we observe the scaled perturbed counterion 
concentration in the direction of the field {9 ~ 0) along 
the cell. This quantity is responsible of the charge po- 
larization contribution as we can see in eqn p.8p . Solid 
black lines, solid colored lines and dashed lines, repre- 
sent the results of the PL, FIS and FIS+L models, re- 
spectively. Different colors stand for different counterion 
sizes. The perturbed counterion concentration is scaled 
as 



(5n*(r) 



eQErsE 

where, using eqns ()2.13|) . (j2.17p and (|2.18p . we have 



(3.14) 



Zee 



n°ir)iY{r) - c/jcir)) Ecos9 (3.15) 



In the inset of Fig. [5] we represent (5n* versus r/a 
in the neighborhood of the particle with the FIS model 
for n™°^ = 15, 18, 22, 25 and 30 M, respectively, from 
bottom to top, just to clarify the transition between the 
PL and FIS cases. 
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100 




1.04 1.06 

r/a 

FIG. 8: Scaled perturbed counterion concentration in the di- 
rection of the field {9 = 0) along the cell for different ion 
sizes. Different lines have the same meaning than those in 
Fig. [2] In the inset, 5n* versus r/a in the neighborhood of 
the particle surface with the FIS model. Colored lines stand 
for 71™"^ = 15, 18, 22, 25 and 30 M, respectively, from bottom 
to top. 



We observe that the perturbed counterion concentra- 
tion decreases asymptotically to zero at the outer surface 
of the cell [h/a = 1.26) in all cases according to Shilov- 
Zharkikh-Borkovskaya boundary condition, eqn (|2.33p . 
In the PL case, we obtain an excess of counterions at the 
rear of the particle, 6* = 0, and a defect of counterions 
at the front of the particle, 9 = it, due to a counterions 
migration from the front to the rear of the particle when 
the external electric field is applied. This excess of coun- 
terions is mainly located very close to the particle surface 
and generates an electric dipole moment that points to 
the direction of the electric field. 

When we take into account the finite counterions size, 
FIS model, the excess of counterions is lower near the 
particle surface in comparison with the PL case, see Fig. 
Ill As we noted in a previous work (Fig. 2 of ref. f40|). 
the equilibrium concentration profile shows a counterions 
condensate that increases its width upon increasing the 
size of the counterions because of steric reasons. There- 
fore, when the external electric field is applied, the excess 
of counterions cannot be located in the condensate, full 
of counterions mainly for high particle charges. This is 
the reason why the main region of excess of counterions is 
now found at farther distances from the particle surface. 
When we consider the excluded region in contact with 
the particle, FIS-I-L model, the perturbed concentration 
profile is shifted to larger distances from the particle sur- 
faces as we see in Fig. [5] In all the cases studied, the 
charge contribution to the induced electric dipole mo- 
ment points to the direction of the electric field which 
penalizes the movement of the particle by the relaxation 
effect (3. 

According to Table HI the charge polarization contri- 



bution is larger the larger the ion size, and even more if 
we take into account the distance of closest approach of 
the counterions to the particle surface. This fact is in 
concordance with an increase in the charge contribution 
to the induced electric dipole moment that we can asso- 
ciate to the charge redistribution depicted in Fig. [S] due 
to the excluded volume effect. 

Summarizing, we have seen that when we introduce the 
ion size effects there is a remarkable diminution of the to- 
tal electric force acting on the particle (see Table H]). This 
decrease grows with the ion size and with the inclusion of 
the excluded region. The diminution of the total electric 
force must be accompanied, in the stationary state, by 
a corresponding diminution of the total hydrodynamic 
force. According to the theory of classical electrokinetics 
3, the total hydrodynamic force could be decomposed 
in a viscous drag and a electrophoretic retardation con- 
tribution. In this frame, the above mentioned decrease of 
the electric body force as ions size effects are considered 
provokes a diminution of the electrophoretic retardation 
force, which opposes to the movement of the particle, 
as it happens in the FIS and FIS-fL cases. However, a 
complete explanation of the observed increase in the elec- 
trophoretic mobility will force us to study the transient 
regime after the application of the external field and the 
evolution of the different forces involved until the station- 
ary state is reached. In the case we are concerned in this 
work, the final result is that the electrophoretic mobility 
of the particle is higher in the FIS and the FIS-I-L cases 
in comparison with the PL model using 11+ ions as coun- 
terions, being these ions commonly found in many exper- 
imental salt-free suspensions. For other ionic species the 
behaviors observed can be different, depending on the 
diffusion coefficient of the ion, because the various con- 
tributions to the total force will be altered. The infiuence 
of the mobility of the ions will be addressed in a future 
work in connection with the corresponding experimental 
results. 



B. Electrical conductivity 

The electrical conductivity of salt-free concentrated 
suspensions with point-like ions behaves classically as fol- 
lows: the conductivity increases as surface charge density 

For large particle 



increases for any volume fraction [21 



charges, the conductivity tends to reach a plateau be- 
cause of the classic counterion condensation effect fio} : 
it appears that a critical particle charge density exists 
beyond which there is no appreciable influence on the 
conductivity. Once this critical charge value is attained, 
increasing the amount of counterions by raising the sur- 
face charge even more simply feeds the condensation re- 
gion, where a high accumulation of counterions takes 
place close to the particle surface, leaving the charge and 
potential outside that region virtually unchanged. There 
is also a conductivity enhancement with particle volume 
fraction because the increasing number of the double- 
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FIG. 9: (a) Electrical conductivity against the surface charge 
density for different particle volume fraction values, (b) Elec- 
trical conductivity against the particle volume fraction for 
different surface charge densities. Solid lines show the results 
for point-like ions. Dashed lines show the results of the FIS+L 
model. 



layer mobile ions is not offset by the presence of the non- 
conducting volume occupied by the particles in the unit 
volume. These behaviors are shown in solid colored lines 
in Fig. inland in solid black lines in Figs. [TUlandITT] 

When we take into account the finite size of the coun- 
terions (FIS model, solid colored lines in Figs. [TU] and 
[TT|) and the distance of closest approach of the counteri- 
ons to the particle surface (FIS+L model, dashed lines in 
Figs. [51 [TUl and [TT|) we observe similar behaviors but the 
numerical values of the electrical conductivity are always 
higher for any counterion size for moderate to high par- 
ticle charges in concentrated suspensions. If the counte- 
rion size approaches to zero, or equivalently n™"^ — oo, 
the results of the FIS and FIS-fL models approximate 
to those of the PL model. Also, for low particle charges 
and low particle volume fractions the results of the three 
models are almost coincident. 



FIG. 10: Electrical conductivity against the particle sur- 
face charge density for different ion sizes, considering (dashed 
lines) or not (solid lines) the excluded region in contact with 
the particle. Black lines show the results for point-like ions. 
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FIG. 11: Electrical conductivity against the particle volume 
fraction for different ion sizes, considering (dashed lines) or 
not (solid lines) the excluded region in contact with the par- 
ticle. Black lines show the results for point-like ions. 



According to Fig. [7] the counterions fluxes not in the 
immediate vicinity of the particle have been enhanced in 
the FIS and FIS+L model in comparison with the PL 
case. This enhancement gives rise to a larger conduc- 
tivity of the suspension, because it predominates over 
the larger PL counterions fluxes very close to the par- 
ticle surface. As these augmented counterions fluxes in 
the FIS and FIS+L models grow with the counterion size 
and with the inclusion of the excluded region, the electric 
conductivity increases as well (see Figs. [9l [TOl and [TT|) . 

Although all the conductivity calculations have been 
performed with 11+ counterions, we have checked that the 
conductivity behavior shown before maintains for other 
counterions species with different diffusion coefficients. 



IV. CONCLUSIONS 

In this work we have studied the influence of finite 
ion size corrections on the electrophoretic mobility of 
spherical particles and the electrical conductivity in salt- 
free concentrated suspensions with 11+ added counteri- 
ons, although the model is also valid for different ionic 
species. The resulting model is based on a mean-field ap- 
proach that has reasonably succeeded in modeling elec- 
trokinetic and rheological properties of concentrated sus- 
pensions. We have used a cell model approach to account 
for particle-particle interactions, and derived a DC elec- 
trokinectic model which include such ion size effects. The 
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theoretical procedure has followed that by Carrique et al. 
[2H with the additional inclusion of finite size counteri- 
ons [iOl and an excluded region of closest approach of the 
ions to the particle surface [37| . 

The results have shown that the finite ion size effect 
has to be taken into account for moderate to high parti- 
cle charges in concentrated suspensions, and even more 
if a distance of closest approach of the ions to the par- 
ticle surface is considered. In the common case of H+ 
counterions, we have found a larger increase of the elec- 
trophoretic mobility and the suspension conductivity the 
larger the ion size. The effect is more noticeable if we 
take into account an excluded region free of ions in con- 
tact with the particle. 

The DC electrokinetic model presented in this pa- 
per will be used to develop a theoretical model on the 
response of a salt-free concentrated suspension to an 
AC electric field including ion size effects. Experimen- 
tal results concerning the DC electrophoretic mobility, 
dynamic electrophoretic mobility, electrical conductivity 
and dielectric response with different counterions species 
should be compared with the predictions of the latter 
models to test them. To perform such comparisons, con- 
centrated suspensions of highly charged particles are re- 



quired, which classically have been very difficult to syn- 
thesize. This is probably the reason that explains the 
lack of experimental studies that could be used to test 
our predictions. Nowadays, existing highly charged sul- 
fonated polystyrene latexes could be good candidates. 
As the authors shown in a previous paper [51[, realis- 



tic salt-free suspensions which include water dissociation 
ions and those generated by atmospheric carbon diox- 
ide contamination, in addition to the added counterions 
released by the particles to the solution, should be con- 
sidered to get closer to the experimental results. The 
inclusion of these chemical reactions in an electrokinetic 
model is not trivial and require a proper electrokinetic 
model that extends the present one. These theoretical 
and experimental tasks will be addressed by the authors 
in the near future. 
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